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ABSTRACT 


The  errors  associated  with  the  analysis  of  acoustic  signals 
generated  by  underwater  explosive  souices  are  analyzed  in 
this  report.  Thus,  the  effects  of  analog  reproduction  and 
subsequent  digital  conversion  of  the  data  on  the  estimation 
of  propagation  loss  are  considered.  In  addition  the  tech- 
niques of  the  spectral  analysis  as  applied  to  these  data 
are  discussed.  It  is  shown  that  the  digitizing  processing 
can  be  repeated  on  a day-to-day  basis  to  the  extent  that 
the  energy  in  frequency  bands  for  either  shots  or  noise 
have  standard  deviations  of  1 to  2$.  The  times  at  which 
the  automatic  shot  processor  detected  a shot  showed  no 
fluctuations  other  than  those  to  be  expected  from  the 
temporal  resolution  imposed  by  the  sampling  rate.  The 
standard  deviations  due  to  repeated  processing  of  trans- 
mission loss  for  a given  frequency  band  is  of  the  order  of 
0.1  to  0.2  dB. 
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I . INTRODUCTION 

This  memorandum  gives  the  results  of  investigations  conducted  to 
obtain  measurements  of  errors  associated  with  various  components  of 
the  processing  system  used  at  Applied  Research  Laboratories,  The 
University  of  Texas  at  Austin,  for  the  analysis  of  ACODAC  analog  data. 
The  significant  components  of  the  analysis  procedure  are  investigated 
individually  to  determine  the  magnitude  of  error  that  each  contributes 
to  the  final  estimate  of  ambient  noise  level  and  propagation  loss.  The 
areas  considered  in  the  following  sections  are  (l)  the  consistency  of 
the  analog- to- digital  conversion  and  shot  processing,  (2)  the  spectral 
analysis  methods  used  to  determine  the  energy  content  of  acoustic 
signals  in  selected  frequency  bands,  (3)  the  effects  of  sampling 
without  phase  control,  and  (4)  the  stability  of  the  ACODAC  calibration 
cignal  as  a function  of  time . The  concluding  remarks  are  made  in 
Section  VI. 


ii.  efhectg  of  tape  playback  /tiro  analog-to- digital 

CONVERSION  UPON  SHOT  PROCESSOR 

To  examine  the  effects  of  analog  tape  playback  and  analog- to- digital 
conversion  of  ACODAC  data  upon  the  shot  processor,  a segment  of  data  was 
subjected  to  repeated  digitization  and  analysis. 

A 1 h segment  of  ACODAC  recording  of  the  CHURCH  ANCHOR  exercise 
was  selected;  the  data  record  contained  approximately  30  shots  plus  one 
1 calibration  signal  sequence.  Topical  envelopes  of  the  600  ft,  300  ft, 
and  60  ft  shots  are  shown  in  Fig.  1.  On  three  different  days,  the  same 
1 h segment  of  recorded  data  was  digitized  five  times  for  a total  of 
15  digitized  records.  These  were  regarded  as  an  ensemble  of  noise 
corrupted  signals;  the  noise  in  this  case  was  due  to  fluctuations 
introduced  by  the  playback  system.  The  digitized  records  were  then 
processed  with  the  ARL  shot  processor,  which  automatically  detects  a 
shot  and  then  does  spectral  analysis  of  the  shot  and  of  a segment  of 
ambient  noise  preceding  it.  The  processor  output  record  includes  the 
energy  of  the  shot  and  the  noise  in  six  frequency  bands  (l/3  octave 
centered  at  12.5;  2p,  30;  100,  l60,  and  200  Hz);  these  are  used  to 
compute  transmission  loss. 

For  each  shot,  the  mean  and  standard  deviation  of  the  computed 
energies  were  estimated;  these  statistics  were  calculated  for  each  day 
(groups  of  5 digitizations)  as  well  as  for  the  ensemble  of  15. 

As  an  example,  consider  the  50  Hz  band  of  the  signals  received  on 
^hydrophone  E from  shots  detonated  at  nominal  depths  of  60  ft.  During 
day  21,  hour  23  (20  December  1972,  IlOO)  of  the  exercise,  the  Bartlett 
was  approximately  500  run  from  station  A;  the  60  ft  shot  arrived 
approx Lmately  20  sec  into  the  13th,  lflth,  23rd,  etc.,  min. 
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FIGURE  1 

REPRESENTATIVE  RECEIVED  SHOT 
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In  Fig.  2 the  energy  (S+N)  averaged  over  the  5 digitizations  per 

day  is  plotted  for  the  digitization  days  versus  minute  of  reception; 

the  energy  scale  is  in  arbitrary  units,  and  the  calibration  signals 

have  not  been  taken  into  account  for  this  figure.  The  largest 

standard  deviation  among  the  digitization  sequences  on  a particular 

day  for  a given  shot  was  0.032  units,  about  the  size  of  the  solid 

dot  symbols.  This  figure  serves  mainly  to  point  up  the  possible 

variation  from  day  to  day  of  the  level  reproduced  from  a given  recording. 

* 

The  shot  detonation  depths  as  provided  by  Underwater  Systems,  Inc., 
are  also  shown  on  Fig.  2. 

Figure  3 shows  daily  averages  for  each  of  the  3 digitization  days 
of  10  Log(S+N-N)  minus  10  Log( Calibration)  versus  time.  Except  for 
source  level  and  hydrophone  sensitivity  compensation,  this  is  effectively 
transmission  loss.  Though  at  most  times  the  differences  are  scarcely 
perceptible,  at  minute  18  there  is  an  almost  0.2  dB  spread.  The 
standard  deviation  of  transmission  loss  for  individual  shots  was 
0.07  to  0.15  dB.  As  may  be  seen  proper  use  of  the  calibration  signal 
removes  the  spread  which  might  be  anticipated  from  Fig.  2.  The 
signal- to-noise  ratios  of  the  signals  used  to  compute  Fig.  3 are  shown 
in  Fig.  k. 

The  300  ft  shot  was  detonated  1 min  before  the  60  ft  shot. 

Figures  5 and  6 show  transmission  loss  curves  similar  to  Fig.  3 for 
the  300  ft  shot  in  the  50  Hz  and  200  Hz  frequency  bands  (absolute  level 
comparison  between  Fig.  3 and  Figs.  5 and  6 should  not  be  made).  The 
small  effect  of  system  variability  upon  transmission  loss  shown  by 
Figs.  3,  5,  and  6 was  typical  of  results  from  all  frequencies  and 
depths.  Also  typical  is  the  variability  of  energy  levels  occurring 
between  closely  spaced  shots. 

* "SUS  Shot  Statistics,"  paper  tape  records  mailed  to  ABL  on 
12  November  1973. 

**N  is  the  noise  energy  level  estimated  from  a portion  of  the  record 
preceding  the  time  of  detection. 
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FIGURE  2 

AVERAGE  SHOT  ENERGY  COMPUTED  FROM 
DATA  DIGITIZED  ON  THREE  DIFFERENT  DAYS 
HYDROPHONE  4,  SIGNAL  PLUS  NOISE 
FREQUENCY  BAND:  1/3  OCTAVE  CENTERED  AT  50  Hz 
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FIGURE  3 

RELATIVE  TRANSMISSION.  LOSS  COMPUTED  FROM 
DATA  DIGITIZED  ON  THREE  DIFFERENT  DAYS 
DEPTH  OF  SHOT  = 60  ft 

FREQUENCY  BAND:  1/3  OCTAVE  CENTERED  AT  50  Hz 

PLOTTED  POINTS  ARE  10  LOG( SIGNAL-NOISE) -10  LOG(CALIBRATION  SIGNAL) 

HYDROPHONE  4 
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FIGURE  4 

AVERAGE  SIGNAL  TO  NOISE  RATIO  COMPUTED 
FROM  DATA  DIGITIZED  ON  THREE  DIFFERENT  DAYS 
SIGNAL  TO  NOISE  RATIO  = 10  LOG[(S+N)/N] 
FREQUENCY  BAND:  1/3  OCTAVE  CENTERED  AT  50  Hz 
HYDROPHONE  4 
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FIGURE  5 

RELATIVE  TRANSMISSION  LOSS  COMPUTED 
FROM  DATA  DIGITIZED  ON  THREE  DIFFERENT  DAYS 
DEPTH  OF  SHOT  = 300  ft 

TRANSMISSION  LOSS  = 10  LOG( SIGNAL-NOISE) -10  LOG (CAL I BRAT ION  SIGNAL) 
FREQUENCY  BAND  = 1/3  OCTAVE  CENTERED  AT  50  Hz 
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FIGURE  6 

RELATIVE  TRANSMISSION  LOSS  COMPUTED  FROM 
DATA  DIGITIZED  ON  THREE  DIFFERENT  DAYS 
DEPTH  OF  SHOT  = 300  ft 

TRANSMISSION  LOSS  = 10  LOG( SIGNAL-NOISE) -10  LOG( CALIBRATION  SIGNAL) 
FREQUENCY  BAND  = 1/3  OCTAVE  CENTERED  AT  200  Hz 
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Examination  of  the  statistics  from  all  hydrophones  and  bands 
showed  the  following. 

1.  On  a single  day,  the  energy  in  the  frequency  bands  for  either 
shots  or  noise  had  standard  deviations  of  1$  to  2$.  The  energy  in  the 
same  bands  of  the  calibration  signals  showed  a similar  spread. 

2.  The  signal  levels  as  digitized  showed  variations  from  day  to 
day  which  yielded  5$  to  10$  (5$=0.2  dB)  deviations  of  energy  in  the 
bands  before  the  calibration  signal  was  used  for  compensation.  This 
day-to-day  variation  was  of  the  nature  of  a daily  change  in  overall 
playback  gain,  though  it  likely  is  due  to  variations  in  the  analog 
tape/playback  head  combination;  the  entire  1 h record,  including 
calibration  signals,  was  similarly  affected,  thus  allowing  compensation 
to  be  made. 

3.  The  times  at  which  the  automatic  shot  processor  detected  a 
shot  showed  no  fluctuations  other  than  those  to  be  expected  from  the 
temporal  resolution  imposed  by  the  sampling  rate.  The  standard 
deviations  of  detection  times  for  individual  shots  were  in  the  range 
5 msec  to  20  msec . The  data  sampling  rate  was  600  Hz,  but  only  every 
fourth  sample  was  considered  for  detection;  from  this,  one  should 
expect  an  rms  jitter  of  13  msec . 

4.  The  standard  deviations  of  transmission  loss  for  a given 
frequency  band  is  of  the  order  of  0.1  to  0.2  dB. 
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III.  EFFECTS  OF  METHODS  OF  DIGITAL  SICNAL  COMBINATION 
ON  SHOT  SPECTRAL  ANALYSIS 


This  section  provides  detailed  descriptions  of  digital  techniques 
used  for  estimating  the  energy  of  the  shot  in  selected  frequency  bands. 
The  current  goal  of  the  analysis  program  is  that  shot  energy  measure- 
ments presented  for  one-third  octave  bands  obtained  by  digital  analysis 
techniques  be  compatible  with  measurements  previously  made  using  analog 
techniques.  Digital  analysis  methods  are  concerned  with  two  practical 
realities . 

1.  The  typical  received  shots  analyzed  in  the  LRAPP  experiment 
can  have  time  durations  of  50  sec  or  1800  samples  at  a sample  rate  of 
600  Hz.  (Durations  up  5 to  4 times  this  length  have  been  observed.) 

A single  Fourier  transformation  of  this  length  is  not  practical  with 
most  computer  capabilities.  This  necessitates  segmenting  the  shot 
record  into  manageable,  contiguous  segments  whose  Fourier  transforms 
are  easily  computed.  This  process,  then,  presents  the  second  problem. 

2.  The  individual  Fourier  transforms  must  be  summed  in  some 
manner  to  give  the  final  power  spectrum  of  the  shot  in  a given 
one-third  octave  frequency  band.  The  summation  can  be  made  using  the 
complex  components,  where  the  phase  information  of  the  Fourier 
transforms  of  the  contiguous  segments  Is  preserved,  or  the  summation 
can  be  made  using  the  energies,  where  the  magnitudes  of  the  individual 
spectra  are  preserved. 

It.  this  section,  a detailed  description  of  these  two  techniques 
is  given.  First  the  theoretical  aspects  of  the  problem  are  considered. 
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A.  Theoretical  Description  of  Digital  Spectral  Analysis 

Consider  a time  function  x(t)*  defined  over  the  interval  (0,T). 
The  Fourier  transform  of  x(t)  is 


'(f)  -<$r|J(t)|  . _£*S(t) 


-icot 

e dt  , 


(1) 


where  o>=2jtf.  Since  x(t)  is  nonzero  only  over  the  interval  (0,T) 


F(f)  - f^xi t) 
J r> 


-hot 

t*  dt 


(2) 


The  spectral  intensity  density  of  x(t)  is 


W(f)  = |F(f)|‘ 


(3) 


The  total  energy  in  the  signal  is 


E 


which  is  also  available  from 


E 


£- 

Tilde  used  to  denote  continuous  representation. 


JP? 


- f x2(t)dt  , 

Jr- 


(*) 


the  spectral  density 


= f W(f  )df 

cr 


(5) 
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The  energy  in  a selected  frequency  band  (f  ,fg)  is 

E12=2y*f2  W(f)df  . (6) 

fl 

The  analog  to  Eq.  (2)  for  discrete  processing  is 


F,  =WKY*  e~(2nlJmA) 
j Z m 
m5) 


where 

x = x(mAt),mAt  e [0,T] 
m 


(7) 


U is  a weighting  factor  to  preserve  Parseval's  relation  (Eqs.  4 and  5). 
It  is  well  known  that 

Fj  * I F<fm> 

m.=-«< 

where  fm=(j+Km)/At.  If  the  sample  rate  is  adequate, 

F(fm)  « 0 , a/0 


sc  that 


F ~ F(j/At) 

Analogous  to  Eq.  6,  the  energy  in  a frequency  interval  [f^fg]  is 
defined  as 
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E12  = 2 


Ir, 


A practical  limitation  to  the  computing  of  Eq.  ^ is  the  size  of 
available  computer  storage.  To  maintain  reasonable  data  cample  length 
it  will  at  times  be  necessary  to  obtain  the  spectral  energy  in  the 
frequency  bands  of  interest  by  transforming  segments  of  the  data  and 
then  accumulating  the  transformation  of  the  segments. 

Let  us  assume  that  the  time  interval  of  the  signal  has  been  broken 
into  N segments  of  length  L;  L=K/N.  We  define 


X ~ X,_i  IT 

m m+iL 


t • 0,1,...,  N-l 


m = 0,1, . . .,  L-l 


and 


L-l 


I -(&tlJm/L) 


Pj' . w y xm‘  , j = o,x,..„  t 

m=t) 


-1 


The  spectral  energy  density  at  frequency  j/At  for  the  signal  in  interval 
I would  be 


V-I'i1'8 


There  are  two  straightforward  accumulations  of  the  frequency  information 
from  the  segments  to  represent  the  spectrum  of  the  entire  signal: 
summation  of  the  complex  components  and  summation  of  the  energies. 


First,  complex,  or  "coherent,"  summation  is  defined  by 


<P 


i 


\ 


It  can  be  shown  that  the  <p  are  exact  samples  of  the  spectra  Fm  for  the 
frequencies  m such  that  m=jN.  For  these  frequencies,  the  effect  in  the 
frequency  domain  of  the  exact  phasing  between  eigenrays  which  occur 
in  different  segments  of  the  record  is  preserved. 

However,  the  spectral  levels  at  the  unsampled  frequencies  m such 
that  jN<m<jN+N  do  not  contribute  to  the  spectrum  as  estimated  via 
coherent  summation.  One  consequence  of  this  is  that  Parseval's 
relation  is  not  preserved.  This  "coherent"  summation  has  an  extra 
advantage  that  the  order  of  summation  and  transformation  may  be 
reversed.  The  energy  in  interval  [F^,F^]  as  estimated  via  coherent 
summation  would  be 


An  alternative  to  the  complex  accumulation  is  summation  of  energies, 
or  "incoherent"  summation.  Under  this  scheme,  the  energy  at  frequency 
3 is  given  by 


17 


It  can  be  shown  that  Parseval's  relation  is  preserved  under  this 
definition.  Unlike  the  complex  summation,  the  spectral  levels  at 

A 

the  unsampled  frequencies,  jN<m<jN+N,  do  contribute  to  E ; the  details 
of  this  interpolation  are  not  easily  described.  However,  the  phasing 
between  eigenrays  in  different  segments  of  the  record  has  no  effect 
upon  . 


B.  Effects  of  Data  Segment  Length  on  Digital  Signal  Summation 


It  is  anticipated  that  the  two  summation  methods  would  begin  to 
differ  with  increasing  magnitude  as  the  number  of  data  segments  is 
increased,  by  decreasing  the  size  of  each  segment  of  a given  detected 
shot  signal.  This  is  illustrated  by  an  example  involving  three 
received  shot  signals.  The  first  two  signals  are  shown  in  Fig.  7 
where  the  upper  trace  corresponds  to  hydrophone  2 and  the  lower  trace 
corresponds  to  hydrophone  3 of  the  ACODAC  system  deployed  at  Site  D 
of  the  CHURCH  GABBRO  propagation  experiment . The  signals  shown  here 
were  received  with  hydrophones  2 and  3 with  the  explosive  source 
206  run  away.  The  ACODAC  time  was  day  9,  hour  20,  minute  14, 
second  30  (8  December  1972,  08l4). 


The  length  of  the  shortest  data  segment  is  512  samples,  which 
corresponds  to  a t't>\  interval  of  0.85  sec  at  a sample  rate  of  600  Hz. 
The  total  interval  ~f  the  shot  analyzed  for  its  power  in  the  l/3  octave 
frequency  band  centered  at  50  Hz  is  32  X 512  = 16384  samples  (27.2  sec). 
In  the  most  extreme  case  a total  of  32  continuous  segments  (0.85  sec 
each)  were  Fourier  transformed  and  summed  coherently  and  incoherently. 

In  the  least  extreme  case  the  entire  16384  samples  (27.2  sec)  were 
transformed  as  one  record.  Between  these  two  extremes  the  number  of 
data  segments  used  were  2,  4,  8,  and  16  (13.6,  6.8,  3.4,  and  1.7  sec 
each  respectively). 


The  reference  for  the  comparison  of  the  coherent  and  incoherent 
summations  is  the  output  of  a digital  recursive  filter  with  Butterworth 
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FIGURE  7 

SHOT  SIGNALS  USED  TO  COMPARE  ENERGY  COMPUTATIONS 


filter  characteristics.  In  the  example  shown  here  a 6- pole  Butterworth 
filter  was  selected  to  match  analog  filters  previously  used  for  the 
analysis  of  ACODAC  data. 

The  frequency  response  of  the  recursive  filter,  shown  in  Fig.  8, 
was  computed  by  generating  sine  waves  of  varying  frequency  and 
measuring  the  output  power.  The  impulse  response  was  also  computed 
and  is  plotted  in  Fig.  9* 

Figure  10  shows  both  the  coherent  and  incoherent  energy  summations 
as  a function  of  the  number  of  segments  transfoimed.  The  coherent 
summation  begins  to  diverge  from  the  reference  energy  as  the  number 
of  segments  increases.  When  the  number  of  data  segments  is  32, 

(each  segment  0.85  sec)  the  percentage  error,  defined  as 

Percentage  Error  = SSS3S.  - .reference  energy  [ 

reference  energy  * 

is  28$.  The  corresponding  decibel  error  is 


Decibel  Error  = 10 


log 


summed  energy 
reference  energy 


-1.4  dB. 


Whenever  four  segments  are  used,  the  corresponding  errors  are  1.6$ 
and  -0.01  dB.  Note  that  the  incoherent  summation  gives  energy  values 
very  close  to  the  reference  energy.  Figure  11  has  the  corresponding 
error  analysis  for  the  received  signal  - ydrophone  3-  For  the 
coherent  summation  case,  the  fluctuations  with  increasing  rms  levels 
will  occur  as  the  size  of  the  segment  is  decreased,  and  these  errors 
will  appear  as  larger  deviations,  either  positive  or  negative,  about 
the  true  power  spectrum.  Thus,  if  other  shot  signals  were  analyzed 
in  this  manner,  it  is  expected  that  the  magnitude  of  the  error  would 
be  similar,  for  a given  segment  size,  but  that  the  error  would  be 
either  positive  or  negative.  This  fluctuation  about  the  true  power 
spectrum  can  be  seen  by  comparing  Figs.  10  and  11.  Note  again  that 
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FIGURE  8 

FREQUENCY  RESPONSE  OF  1/3  OCTAVE 
RECURSIVE  FILTER 
(6 -POLE  BUTTERWORTH) 
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FIGURE  10 

COMPARISON  OF  1/3  OCTAVE  BAND  ENERGIES 

OBTAINED  BY  COHERENT  AND  INCOHERENT  SUMMATION 

SHOT  1 : CHURCH  GABBRO  HYDROPHONE  2 

FREQUENCY  RANGE:  44.5-56.1  Hz 
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-•COHERENT  SUMMATION 
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SIGNAL-TO-NOISE  RATIO  OF  FILTERED  SH0T«7  dB 


AS-74-375 


NUMBER  OF  SEGMENTS 
4 8 


li  |i  1 1 ,11.1 ! I 


Ilii  i 

$ iw 

il 

TfrTTTjT 

I » • [ t*  ♦ 

Jjj?  | 

it# 

jijj  ; 

ii  In  i hi 

* ■ * i ♦ > i'i‘ 

7*  ■W'-'rrt 
f I * ! 

ill; ! 

i ' j \ j 

4r  -ri:m44 


firt'F  /t'ii'Hwi 

/13.6 


ijfi 

nil 

IhT 

my 

tnr 

l 

1 

i 
♦ * ! 1 

jlil 

.[]] 

:iii 

jilj 

lilt 

Ul: 

Hi: 

Sffi 

m;  7j 

jfljtj  1 i 

; u . 

till 

;u; 

i!ii 

i i • * 

mi  ‘ i 

14 

Ilf! 

iiilGl]] 

it  j i : 
iii:  N 

pn 

i ■ 1 1 . i j 

i :ii:  1!  : 
:t  rrtr  Hr; 

1 ;ll!  Hi 

1 I l . 1 1 ! 1 . i 

[i 

:|H 

iUU 

• i i i . • 

in: 

i » ! j i ♦ • t 

» l » » • * • u 
: 

Mj!  j|i!  i li 


|1  |lli 

t!:  t!H  tiH  | i:  :•  ! t 


i « • i m ♦ • i » ' ] i rnr 

nHMPlij 


13.6  6.8  3.4  1 .7  S .85 

FFT  1ength  ^sec)  PERCENTAGE  ERROR  =41% 

(Total  data  length  27.2  sec)  dB  ERROR  = -2.3  dB 


PERCENTAGE  ERROR  =9.8% 
dB  ERROR  = 0.41  dB 


FIGURE  11 

COMPARISON  OF  1/3  OCTAVE  BAND  ENERGIES  OBTAINED 

BY  COHERENT  AND  INCOHERENT  SUMMATION 

SHOT  1 : CHURCH  GABBRO  HYDROPHONE  3 
FREQUENCY  RANGE:  44.5-56.1  Hz 
— REFERENCE  LEVEL 

...  COHERENT  SUMMATION 
000  INCOHERENT  SUMMATION 

SIGNAL -TO-NOISE  RATIO  OF  FILTERED  SH0T^7  dB 
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incoherent  summation  gives  results  essentially  equivalent  to  the 
reference  level  for  any  number  of  segments  up  to  32.  From  these 
results  it  may  be  concluded  that,  for  data  of  this  type,  the  loss 
of  phase  coherence  (incoherent  summation)  between  record  segments 
causes  less  error  than  does  the  loss  of  information  at  unsampled 
frequencies  (coherent  summation).  Tho  agreement  between  the  incoherent 
summation  and  the  reference  suggests  that  there  is  no  critical  phasing 
among  components  of  the  signal  separated  by  more  than  a few  seconds . 

One  further  point  was  made  with  these  two  shot  signals  (Fig.  7). 

In  a practical  situation  the  received  shot  signal  is  never  known 
exactly  due  to  the  influence  of  background  noise  (ambient  or  system). 

A question  arises  concerning  the  effect  that  an  ambiguity  in  shot 
location  has  on  the  estimation  of  the  energy  in  a given  frequency  band 
with  coherent  or  incoherent  summation.  Some  insight  is  obtained  by 
considering  Fig.  12,  which  is  identical  to  Fig.  ^ except  that  the 
beginning  of  the  defined  shot  is  earlier  by  a time  of  approximately 
3.5  sec.  Analyzing  this  signal  in  precisely  the  same  manner  as  the 
previous  example  (Figs.  10  and  11 ) the  data  in  Figs.  13  and  Ik  were 
generated.  As  can  be  seen  the  effect  of  not  knowing  the  epoch  exactly 
resulted  in  approximately  0.1  dB  difference  in  the  estimation  of  the 
shot  energy  by  coherent  or  incoherent  summation.  It  is  expected, 
however,  that  the  ambiguity  in  epoch  will  have  more  serious  consequences 
as  the  signal- to-noise  ratio  of  the  shot  decreases  to  a relatively 
lower  level. 

The  recorded  shot  signals  shown  in  Fig.  7 ho  not  exhibit  any 
strong  multipath  structure  which  conceivably  could  yield  different 
results  with  respect  to  the  comparisons  of  coherent  and  incoherent 
summations.  A signal  with  a strong  multipath  structure  evident  is 
shown  in  Fig.  15.  This  signal  was  received  with  hydrophone  2 of  the 
ACODAC  system  located  at  Site  D with  a range  of  104  nm  from  the  source. 
The  ACODAC  time  was  day  10,  hour  7,  minute  23,  second  10.  A comparison 
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SHOT  SIGNALS  USED  TO  COMPARE  ENERGY  COMPUTATIONS 
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FIGURE  13 

COMPARISON  OF  1/3  OCTAVE  BAND  ENERGIES  OBTAINED 
BY  COHERENT  AND  INCOHERENT  SUMMATION 
SHOT  2:  CHURCH  GABBRO  HYDROPHONE  2 
FREQUENCY  RANGE:  44.5-56.1  Hz 
- — REFERENCE  LEVEL 
• • • COHERENT  SUMMATION 
000  INCOHERENT  SUMMATION 
SIGNAL-TO-NOISE  RATIO  OF  FILTERED  SH0T=*7  dB 
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dB  ERROR  = 0-51  dB  dB  ERROR  = -2.1  dB 


FIGURE  14 

COMPARISON  OF  1/3  OCTAVE  BAND  ENERGIES  OBTAINED 
BY  COHERENT  AND  INCOHERENT  SUMMATION 
SHOT  2:  CHURCH  GABBRO  HYDROPHONE  3 
FREQUENCY  RANGE:  44.5-56.1  Hz 
— REFERENCE  LEVEL 
. • • COHERENT  SUMMATION 
000  INCOHERENT  SUMMATION 
SIGNAL -TO-NOISE  RATIO  OF  FILTERED  SH0T»7  dB 
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of  the  coherent  and  incoherent  energies  in  the  one-third  octave  band 
at  50  Hz  is  shown  in  Fig.  16,  and  the  errors  are  of  the  same  magnitude 
as  those  of  the  shots  with  the  fused  multipath  structure. 

C.  Comparison  of  Methods  of  Digital  Signal  Summation 

For  the  first  ACODAC  shot  analysis  performed  at  ARL,  both  the 
coherent  and  incoherent  summation  results  were  included.  These 
preliminary  data  from  approximately  700  shots  allow  a statistical 
comparison  of  the  two  techniques . 

The  first  experimental  data  are  summarized  in  Tables  I and  II 
(ACODAC  hydrophones  4 and  5*  respectively)  where  the  means  and 
variances  of  signal- plus- noise  (S+N)  and  signal- to- noise  ratios 
(S+N)/N  corresponding  to  several  one- third  octave  frequency  bands  are 
listed  for  both  coherent  and  incoherent  summation  methods.  Also  listed 
are  the  correlation  coefficients  (covariance)  computed  between  the 
coherent  and  incoherent  forms  of  these  estimates . 

To  estimate  the  shot  energy  S,  which  one  must  do  to  compute 
propagation  loss,  the  number  (S+N)  is  first  computed  and  from  this  is 
subtracted  the  noise  estimate  N.  Similarly  the  signal- to- noise  ratio 
determines  the  accuracy  of  the  estimates  of  important  characteristics 
of  the  shot. 

As  an  example  of  how  the  values  of  Tables  I and  II  are  computed, 
the  mean  of  the  coherent  estimate  of  (S+N)  (Table  II ) is 

722 

" 752  I (S+H)Uh 

i=l 

where  (S+N)£  ^ is  the  estimated  energy  in  a one-third  octave  band  at 
the  center  frequency  f.  The  sura  is  over  722  received  shot  signals. 

The  shots  were  distributed  with  respect  to  the  number  of  4096  sample 


FIGURE  16 

COMPARISON  OF  1/3  OCTAVE  BAND  ENERGIES 

OBTAINED  BY  COHERENT  AND  INCOHERENT  SUMMATION 

THE  SHOT  ANALYZED  IS  SHOWN  IN  FIGURE  15 

FREOUENCY  RANGE:  50  Hz 
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transforms  as:  2 transforms,  38$;  4 transforms,  30$;  5 transforms,  21$; 

6 transfoims,  11$.  Within  a given  shot  the  length  of  a basic  time 
segment  is  4096  samples  (approximately  6.83  sec  at  600  Hz  sample  rate). 
It  is  these  contiguous  time  segments  that  are  Fourier  transformed  and 
added  coherently  to  obtain  the  spectrum  of  the  one- third  octave 
frequency  band.  These  numbers  from  722  shots  are  then  averaged  to 
give  the  mean  of  (S+N)  according  to  the  above  expression. 

The  variance,  likewise,  is 


VAR( S+N ) 


f _1_ 
coh  = 722 


(S+N)f;C0h 


The  correlation  coefficient,  computed  to  determine  to  what  extent  the 
coherent  and  incoherent  sums  are  equivalent,  is  defined  as  (for  the 
(S+N)  estimate) 


1 

722 


C0V(S+N)  = 


722 


i=J 


<S+N>l,coh  -Mcoh 


(S+N)f 

X> 


, . -M  ^ 

incoh  incoh 


VARcoh  (S+K)  VARcoh  (S+N) 


W 


where  and  are  the  corresponding  means  of  (S+N),  coherent 

and  incoherent.  For  the  ( S+N ) /N  estimates,  the  covariance  was  computed 
with  data  expressed  in  decibels.  As  can  be  seen,  the  covariance  is 
normalized  by  the  variances  of  the  estimates.  Thus,  the  maximum 
covariance  is  unity,  a value  which  occurs  if  the  coherent  and  incoherent 
techniques  are  identical. 


As  Tables  I and  II  illustrate,  the  coherent  and  incoherent  summation 
techniques  are  essentially  identical  in  that  most  of  the  covariances 
are  very  close  to  unity,  the  majority  being  greater  than  0.99*  The 
lowest  value  is  0.92^9,  which  is  the  covariance  between  the  coherent 
and  incoherent  estimates  of  (S+N)/N  in  the  one- third  octave  frequency 
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STATISTICS  COMPARING  COHERENT  AND  INCOHERENT  TECHNIQUES  FOR  ESTIMATING 
POWER  SPECTRA  IN  VARIOUS  l/3  OCTAVE  FREQUENCY  BANDS 
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interval  centered  at  12.5  Hz.  Further  evidence  of  the  similarity  of 
the  coherent  and  incoherent  summations  is  given  by  the  relative 
closeness  of  the  corresponding  means  of  the  estimates  of  (S+N)  and 
(S+N)/N  and  the  similarities  seen  in  the  corresponding  variances. 

In  connection  with  this  study  a large  quantity  of  shot 
data  was  analyzed  to  obtain  histograms  of  the  quantities 
10  loeU3+N)coh/(S+N)lncob]  and  10  log[(S+N-a)ooh/(S+N4)lncolj], 
where  S+N-N  is  essentially  the  propagation  loss  of  the  explosive  signal. 

The  number  N is  an  estimate  of  the  background  noise.  These  histograms 
revealed  that 

1.  the  ensemble  means  of  the  quantities  were  mostly  0,1  dB  or  § 

less;  usually  0.01  dB  for  (S+N),  0.1  for  (S+N-N).  That  is  to  say 

( )coh  is  an  unbiased  estimate  of  ( )incoh.  j 

2,  The  rms  difference  varied  with  frequency  and  integration 
time  and  were  usually  less  than  1 dB. 
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IV.  EFFECTS  CF  SAMPLING  WITHOUT  PHASE  CONTROL 


Mechanical  variations  (in  time)  in  both  the  recording  and  playback 
tape  recorders  used  with  the  analog  data  can  result  in  frequency 
deviations  or  an  increase  in  bandwidth  of  the  original  acoustic  data. 
These  variations  can  be  compensated  for  in  part  by  the  recording, 
simultaneously  with  the  acoustic  data,  of  a reference  frequency  track 
which  is  used  on  playback  as  a clock  for  the  analog- to- digital  conversion 
of  the  data.  To  demonstrate  the  effect  of  not  compensating  for  the 
mechanical  variations  in  the  record- playback  process,  the  power  spectra 
of  Figs . 17  and  18  are  shown , Figure  17  shows  the  estimated  power 
spectrum  of  the  100  Hz  ACODAC  calibration  signal  which  was  digitized 
with  the  50  Hz  signal  of  the  time  code  channel  as  a reference  clock 
(multiplied  to  600  Hz).  In  contrast,  Fig.  18  shows  the  corresponding 
power  spectrum  that  is  estimated  from  the  calibration  signal  digitized 
without  the  benefit  of  the  50  Hz  reference  track.  The  spread  in  the 
spectrum  due  to  the  lack  of  phase  control  changes  from  approximately 
0.1  Hz  to  approximately  0.5  Hz  at  the  half  power  point,  illustrating 
the  requirement  of  a reference  track  to  compensate  for  the  mechanical 
variations  in  the  tape  recording- playback  process. 
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FIGURE  17 

POWER  SPECTRUM  OF  100  Hz  CALIBRATION  SIGNAL 
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V.  STABILITY  OF  CALIBRATION  SIGNAL  AS  A FUNCTION  OF  TIME 


The  calibration  signals  are  an  integral  part  of  the  ACODAC  system 
and  their  stability  during  the  course  of  a recording  period  is  a 
prerequisite  for  the  accurate  estimation  of  the  received  acoustical 
energy.  For  this  reason  the  calibration  signals  of  a CHURCH  GABBRO 
tape  were  checked  for  energy  stability  over  a 24  h period.  To  accomplish 
this,  five  3.6  min  calibration  segments,  with  a separation  of  6 h 
between  two  successive  segments,  were  digitized  at  a sample  rate  of 
6 00  Hz.  This  was  done  for  both  the  50  Hz  and  200  Hz  calibration 
signals  on  hydrophones  4,  5>  and  6.  In  addition,  this  was  done  on 
two  SUS  runs . 

Once  in  digital  form,  the  energies  of  each  3.6  min  data  segment 
in  one-third  octave  bands  about  the  center  frequencies  (50  and  200  Hz) 
were  computed.  The  results  are  summarized  in  Table  III,  which  shows 
relative  powers,  expressed  in  decibels.  The  results  shown  here 
indicate  that  the  calibration  signal  is  very  stable;  a typical 
variation  in  power  over  a 24  h period  is  of  the  order  1 dB. 
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VI . SUMMARY 

The  data  presented  in  this  report  describe  the  errors  associated 
with  the  processing  of  the  acoustical  data  obtained  from  the  ACODAC 
receiving  arrays  and  allow  comparison  of  these  errors  with  the 
fluctuations  due  to  the  variations  of  the  acoustical  environment 
(see  Figs.  2 through  6).  Specifically , it  has  been  determined  that 

1.  on  a single  AR1  digitizing  day,  for  repeated  digitization, 
the  energy  in  frequency  bands  for  either  shots  or  noise  had  standard 
deviations  of  1 to  2$;  the  energy  in  the  same  bands  of  calibration 
signals  showed  a similar  spread. 

2.  the  signal  levels  as  digitized  showed  variations  from  day 
to  day  which  yielded  5 to  10$  (5$3D.2  dB)  deviations  of  energy  in  the 
bands  before  the  calibration  signal  was  used  for  compensation; 

J.  the  times  at  which  the  automatic  shot  processor  detected  a 
shot  showed  no  fluctuations  with  repeated  processing  other  than  those 
to  be  expected  from  the  temporal  resolution  imposed  by  the  sampling 
rate  (i.e.,  approximately  20  msec); 

1*.  the  standard  deviations  due  to  repeated  processing  of 
transmission  loss  for  a given  frequency  band  is  of  the  order  of 
0.1  to  0.2  dB. 

5.  the  ensemble  means  of  the  quantity  10  log[(S+N-N)(.  ;oiy(S+N-N)incoh] 
obtained  from  analysis  of  approximately  700  shots  were  of  the  order  of 

0.1  dB  or  less  with  a rms  value  of  less  than  1.0  dB;  these  numbers 
indicate  the  closeness  of  measuring  transmission  loss  by  either  technique. 

6.  to  prevent  an  introduction  of  frequency  smearing  due  to 
mechanical  variations  in  the  recording  and  playback  operations  of  the  tape 
recorder,  a reference  track  should  be  used  in  any  processing  of  the  data, 
whether  analog  or  digital  (see  Figs.  17  and  18). 
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